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Abstract 

Relativistic corrections to the positronium decay rate are calculated. They are 
close to 40(a/7r) 2 and 46(a/7r) 2 for singlet and triplet states respectively. 
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1. The strong disagreement between the experimental value of the orthopositronium 
decay rate[jl| 

K^ S = 7-0482 ± 16 lis' 1 . (1) 
and its theoretical value which includes the order a and a 2 log(l/o;) corrections [@, || |J f| 

= 7.03830 fis' 1 . (2) 

is a serious challenge to modern QED. For the disagreement to be resolved within the QED 
framework the correction ~ (a /it) 2 , which has not been calculated up to now, should enter 
the theoretical result (^) with the numerical factor 250 ± 40 which may look unreasonably 
large. 

Such a hope is not as unreasonable however. The argument is as follows|6[]. The 
large, —10.28, factor at the a/n correction to the decay rate (see (|[)) means that the 
typical magnitude of the factor at the a/ir correction to the decay amplitude is roughly 
5. Correspondingly, this correction squared contributes about 25(a/n) 2 to the decay rate. 
Indeed, numerical calculations |7J have given factor 28.8±0.2 at (a/ir) 2 in this contribution. 

Moreover, it is only natural to expect that the interference of the second-order radiative 
correction to the amplitude with the zeroth-order amplitude should contribute about twice 
as much to the decay rate as the square of the first-order correction. In other words, 
the natural scale for the total second-order radiative correction to the decay rate can be 
about |H 

100(a/7r) 2 . (3) 



-<o— Ps 
th 



ma 6 — - 1-10.28 a 2 log- 

9n 7r 3 a 



A similar conclusion is made in a recent paper starting from the Pade approximants. 

One more class of large contributions to the positronium decay rate is relativistic correc- 
tions. A simple argument in their favour is that the corresponding parameter (v/c) 2 ~ a 2 
is not suppressed, as distinct from that of usual second-order radiative corrections, (a/n) 2 , 
by the small factor 1/tt 2 ~ 1/10. In this article we present the results of calculations of 
relativistic corrections to the positronium decay rate. 



This problem was addressed previously in Refs. [pH, [11]. We differ essentially from those 
authors in the approach to the problem and, which is more essential, in the conclusions 
made. The origin of the disagreements will be elucidated below. 

As to the relativistic correction to the parapositronium decay rate, also obtained in 
the paper, its calculation was started by us as a warm-up exercise for the much more 
complicated orthopositronium problem. However, the correction in the singlet case also 



turns out large, quite close to the sensitivity of the recent experiment [12 



2. The central point when treating the relativistic corrections to the positronium decay 
rate is as follows. Calculating the decay amplitude we have to integrate the annihilation 
kernel over the distribution of the electron and positron three-momenta p. To lowest 
approximation in v/c the kernel, both for para- and orthopositronium, is independent of 
those momenta and we are left with the integral over p of the nonrelativistic wave function 
in the momentum representation which is equivalent to ip(r = 0) in the coordinate one. 
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However, already to first order in (p/m) 2 the momentum integral 



/8 yTFfl^ 
dp (p/m) 2 2 - — (4) 

linearly diverges at p —>■ oo (a = 2/ ma is the positronium Bohr radius). Crucial for the 
problem is the following observation. The true relativistic expression for the annihilation 
kernel does not grow up at p — > oo, as distinct from its expansion in p/m. So, its integral 
with ip(p) in fact converges. Let us transform therefore the integral over \p\ into that from 
— oo to +00 and shift the integration contour into the upper halfplane. We will first come 
across the wave function pole at p — ima/2 and then the relativistic branching point at 
p = im which originates from the amplitude and is not related by itself to the wave function. 
It is obviously just the pole contribution which corresponds to the relativistic correction 
we are looking for. This contribution can be easily calculated and constitutes 

-\c?^r = 0) 

. In other words, the recipe for treating the relativistic corrections originating from the 
amplitude is just to make the substitution 

(p/m) 2 = v 2 ^ -\a 2 (5) 

in them. One may wonder about the sign in rhs of this relation. Let us have in mind 
however that the main contribution to the integral comes from the relativistic cut. That 
contribution corresponds to usual radiative corrections ~ a /it and is of course much larger 
than the effect ~ a 2 we are interested in. 

The next point essential for our consideration is the use of noncovariant perturbation 



theory (see, e.g., Ref.[|T^]) which allows us to treat in a natural way the positronium binding 
energy. 

Let us start with a more simple case of parapositronium. Here the noncovariant anni- 
hilation amplitude can be written as 

M = 4naV+(e 2 a) A ^- k) - k) &5)U + (1 ~ 2) , 



V 




{-^^i 1 -^^^ - )^= x ^^[l + - 7 ^)[ £1.(6) 



In this expression \ and <fi are nonrelativistic spinors; E = 2m — ma 2 / 4 is the positronium 
total energy; ei j2 and &i j2 are the polarizations and momenta of the photons; uj\ = u> 2 = 
u = E/2 are their frequencies; e(p) = ^m 2 + p 2 ; 

1 / ap + 0m \ 

are the projectors onto the positive and negative energy states of a fermion with a momen- 
tum p correspondingly. The Coulomb interaction in the intermediate state can be neglected 
since the momentum of one particle in it is close to m. 
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The expansion of the amplitude in p/m is straightforward. Averaging over the directions 
of p (an 5-state is under discussion) and using relation (H) we obtain 



M + 5M = [l + a 2 (l + y^)}M (7) 
2 o 

where M is the lowest order annihilation amplitude. The corresponding relative correction 
to the decay rate is 

^ = a 2 (l + ^) = 1.35a 2 . (8) 

3. The calculation of relativistic corrections for the triplet state decaying into three 
photons is much more tedious problem. We believe that have managed to simplify it 
considerably, but still it is too lengthy to be presented in detail in a short letter. So, only 
its brief outline is given below. 

The construction of noncovariant perturbation theory amplitude is straightforward. But 
it is convenient to treat it in a slightly different way than it was done for the singlet state. 
First of all rewrite the initial energy E in the perturbative denominators as 

E = E- 2e(p) + 2e(p) 

and expand the amplitude in 

E - 2e(p) 

(here we use recipe (||) for p 2 /m 2 ). 



ma 2 



Zeroth term of the expansion transforms into usual covariant Feynman amplitude for 
electron and positron with 4-momenta (e(p), ±p). Now we expand this "covariant" ampli- 
tude in p/m, average the terms of second order in p/m over the directions of p and make 
the substitution ([5|). The interference of this a 2 -correction with the lowest order amplitude 
after the summation over the photons polarizations and integration over the final phase 
space generates the following correction to the decay rate: 

STc _ 31V-240 

r a 16(vr 2 -9)- 1 } 

This correction is conveniently combined with that originating from the phase space 
correction. It can be easily demonstrated that the shift of the total energy from 2m to 
E = 2m — ma 2 /A changes the phase space and therefore the decay rate by 

= -« 2 -- (10) 



In this way we come to the following total "covariant" correction to the decay rate 

5T C + 5T P 2 27vr 2 - 204 



or — r~n r- ( n ) 

16(tt 2 -9) V ; 
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We have checked that if one goes over from a 2 to v 2 in Eq.(^) according to prescription 
(||) and changes the phase space shifting 2m — > 2m + mv 2 , the result coincides with the v 2 
correction to the probability of 37 annihilation of free electron and positron in 3 Si state 
calculated in Ref. |T0| (see also Ref. [ 



The correction to the decay rate induced by the term of first order in E — 2e(p) = ma 2 /2 
in the expansion of exact amplitude demands numerical calculations which give 

^ = 0.807a 2 . (12) 

Let us note here that the weird term with a/2 in the correction to the singlet decay rate 
(see Eqs. (|7|), (§)) is of the same "noncovariant" origin. 

Going back to the triplet decay rate, we have to note that our results for the a 2 cor- 



rections themselves differ completely from those of Refs. |10|, [TT[| . It is not so much due 
to the "noncovariant" correction ( |T2| ) completely lost there, this correction is not so large 
numerically. The main problem is that of translating v 2 into a 2 . In Ref.||10|| they use the 
prescription v 2 — > a 2 which leads to a wrong sign of the a 2 correction (though to a rea- 
sonable absolute value). On the other hand, the prescription v 2 — > —a 2 /4 used in fact in 



Ref.[flT| gives a correct sign, but strongly underestimate the effect. Quite possibly however, 



there is no direct contradiction between their result and ours, since as it is stated explicitly 



in Ref. , their consideration refers to a part of relativistic corrections only. 



4. Let us consider at last the effects originating from relativistic corrections to the wave 
function ip(r) itself. We will use here the Breit equation following to some extent Ref. ||. 
Para- and orthopositronium can be treated here in parallel. 

The part of the Breit Hamiltonian (BH) that corresponds to the relativistic corrections 
to the dispersion law of the particles and to their Coulomb interaction, 

V e = -/^ + — 2 5(r) , (13) 



a 3 



can be easily transformed to 

v. = & ("> 

Here and below we omit constant terms in the perturbations (obviously, they do not change 
the wave function) and substitute —ma/ 2 for d r acting on the ground state positronium 
wave function. 



The next spin-independent term in the BH 

a 
2m 2 



V m = -77^5- (p 2 + \r(rp)fi) , (15) 



describes the magnetic electron-positron interaction due to the orbital motion. For the 
ground state it transforms into 



The last term in BH of interest for our problem is the contact spin-spin interaction 

V ss = — 2 A5{r); A= 7 -S{S +\) -2. (17) 



It is conveniently rewritten as 
V ss = A 



1 \ a] . a 2 
E 



Am 



r 



+ A -; E = p 2 /m- a/r. (18) 

4mr 2 



The terms a 3 /8r, a 3 /4r from Eqs. fll4|), ( 16|) taken together shift obviously the 
coupling constant a —>■ a(l — 3a 2 /8) which leads to the following relative correction both 
to the | "^(0)1 2 and decay rate 

8T-, 9a 2 

-jr = (19) 

As easily one can calculate the relative correction due to the commutator term in Eq. 
(@): 

Let us turn now to the singular part of the Breit perturbation 

V2 = A; x=a2 (j~b- ( 21 ) 

mr z 4 2 

The normalized solution of the radial wave equation 

fi^r-4 + — + m^i2 = (22) 

is 

R = 2(ma/2) 3/2 [l - A(3 - C)]{mar) x exp [-(1 - A)mar/2] (23) 

where C = 0.577 is the Euler constant. The eigenvalue E deviation from —ma 2 /A is by 
itself irrelevant to our problem. The corresponding relative correction to the |t/>(0)| 2 and 
decay rate constitutes obviously 

5 Jl = _2A[(3 - C) - log(maro)] (24) 

where r ~ 1/m is the distance at which the annihilation takes place. The logarithmically 
enhanced part of this correction 

« 2 log(l/«) | S S Z I ■ (25) 

has been calculated previously for triplet (see formula (0) and singlet cases in Refs. j|] and 
|| respectively. So, we omit it and in this way come to the following total relativistic 
correction due to the ^-function modification: 

ST ±= 2 J 31/8-2C-21og(mr ), 5 = 
T \ -19/24 +1/3C + 1/3 log(mr ), S = 1 ' 1 ' 
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We believe that ±1 is a fair estimate for the scatter of possible values of log(mr ) 
introduced by the uncertainty in the short-distance cut-off r . On the other hand, the 
cut-off of the logarithmic contribution at the atomic distances has been taken care exactly 
in our consideration. 

Our result for the atomic relativistic correction in orthopositronium, (— 19/24+l/3C)a 2 = 
—0.6a 2 , differs from that given in Ref.[ll[]. We cannot explain the disagreement, since the 
authors of Ref. [O present only their numerical result for this correction: 1.16a 2 . 



5. To summarize, the total relativistic corrections in para- and orthopositronium con- 
stitute, respectively, 

= 4.1a 2 = 40(-) 2 , S = 0; (27) 

r 7r 

^ = 4.7a 2 = 46(-) 2 , S = l (28) 

r 7r 

(it is instructive perhaps to present these relativistic corrections in usual "radiative" units 
(a/ir) as well). The terms with log(mr ), omitted here, introduce the errors which we 
estimate as 

±2a 2 = ±20(-) 2 , 5 = 0; (29) 



7T 



±ia 2 = ±3(-) 2 , S = l. (30) 

3 7T 



As to the orthopositronium decay rate, our correction ( f28"D and that of Ref. [[!]], taken to- 
gether, reduce essentially the gap between the theory and experiment, from (250±40)(a/7r) 2 
to (175±40)(a/7r) 2 . 

In parapositronium the magnitude of the calculated correction is close to the present 
experimental accuracy Here there are no special reasons to expect that true radiative 
corrections are as large. So the measurement of the effect looks sufficiently realistic. 
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